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Abstract 

We investigate the possibility to detect the scalar Higgs boson decay H ^ bb m. the 
associated Z and bb production at the LHC using the /cr-factorization QCD approach. 
Our consideration is based on the off-shell (i.e. depending on the transverse momenta of 
initial quarks and gluons) production amplitudes of q*q* —?• ZH —)■ Zq'q', q*q* —)■ Zq'q' 
and g*g* —)■ Zq'q' partonic subprocesses supplemented with the Catani-Ciafoloni-Fiorani- 
Marchesini (CCFM) dynamics of parton densities in a proton. We argue that the H ^ bb 
signal could be observed at large transverse momenta near Higgs boson peak despite the 
overwhelming QCD background, and point out an important role of angular correlations 
between the produced Z boson and 6-quarks. 
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In 2012 , during the search performed by the CMS and ATLAS Collaborations at the 
LHC, the scalar Higgs boson H with a mass mu near 125 GeV has been discovered BE], 
giving us the conhdence in the physical picture of fundamental interactions which follows 
from the Standard Model Lagrangian. Some time later, the ATLAS Collaboration has 
reported hrst measurements of the Higgs boson differential cross sections in the 77 decay 
mode [3]. The measured cross sections were found a bit higher than the central next- 
to-next-to-leading order (NNLO) expectations [IHn] and those matched with soft-gluon 
resummation carried out up to next-to-next-to-leading logarithmic accuracy (NNLL) [TUI 
although no signihcant deviations from theoretical predictions are observed within 
the uncertainties |3]. The signihcant signal was detected also in channels where the Higgs 
boson decays into the ZZ 01 WW pairs nans]. The interaction of the Higgs particle 
with the massive Z and W bosons indicates that, as it was expected, it plays a role in 
electroweak symmetry breaking. However, the interaction with the fermions and whether 
the Higgs held serves as the source of mass generation in the fermion sector still remains 
to be established. Since Higgs boson with mass mriH ~ 125 GeV decays mainly into a 
beauty quark-antiquark pair HI, the observation and study of the H ^ bb decay (which 
involves the direct coupling of the Higgs boson to beauty quarks) is therefore essential in 
determining the nature of the newly discovered boson. 

The most sensitive channel for the H ^ bb events at the LHG is the production of 
Higgs particle in association with the Z boson na. Despite the largest branching fraction 
(~ 58%), the H ^ bb hnal state is a more difficult for the experimental observation 
compared to the signatures provided by the diphoton or diboson decay modes due to 
small signal over background ratio. One of main backgrounds for the associated Higgs 
and Z boson production is the associated production of Z boson and two 6 -quark jets. The 
corresponding cross sections, calculated at the NNLO level (see m), are several orders 
of magnitude larger than the Higgs boson signal. However, recently GMS Gollaboration 
reported [I 6 ] an excess of events above the expected background with a local signihcance 
of 2.1 standard deviations, which is compatible with a Higgs boson mass of 125 GeV. 
Earlier, the GDF and DO Gollaborations at the Tevatron also reported [TU] evidence for 
an excess of events in the 115 — 140 GeV mass range, consistent with the mass of the 
Higgs boson observed at the LHG. 

The experimental searches [T51IT6] are stimulated us to investigate the associated Higgs 
(decaying into a bb pair) and Z boson production as well as corresponding main back¬ 
ground process, associated production of Z boson and two 6 -quark jetfl using the kx- 
factorization approach of QGD [T3|TU]. A detailed description of this formalism can be 
found, for example, in reviews [19]. We only mention that the main part of higher-order 
QGD corrections (namely, NLO -|- NNLO -|- N^LO -|- ... contributions which correspond 
to the logl/x enhanced terms in perturbative series) is effectively taken into account in 
the fc-r-factorization approach already at leading order, and it provides solid theoretical 
ground for the effects of initial parton radiation and transverse momenta of initial quarks 
and gluons. Recently, the fcr-factorization QGD approach was successfully applied [201121] 
to describe the ATLAS data [3] on the inclusive Higgs production in the diphoton decay 
mod^. 

Let us start from a short review of calculation steps. Our consideration is based on the 

^Other background processes, like as tt pair, diboson or QCD multijet production are out of our 
present consideration. 

^In our opinion, the results m suffer from the problem of double counting and contain the wrong 
numerical factor. 
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off-shell (depending on the transverse momenta of initial partons) partonic subprocesses: 


+ q*{k2) ^Z + H^ Zip) + g'(pi) + g'(p 2 ), ( 1 ) 

q*{ki) + q*ik 2 ) ^ Z{p) + g'(pi) + q'{p 2 ), (2) 

9*iki) + g*ik2) Zip) + g'(pi) + q'ip2), (3) 

where the four-momenta of all corresponding particles are given in the parentheses (see 
Fig. 1). The subprocesses (2) and (3) correspond to the main QCD background to the 
associated Higgs and Z boson production. Note, to calculate the production amplitudes, 
we apply the reggeized parton approach [2^123] . which is based on the effective action 
formalism 1211 , currently explored at next-to-leading order I 2 SI, and take into account the 
virtualities of both initial quarks and gluons. In this point our consideration differs from 
the one based on the collinear QCD factorization, where these virtualities are not taken 
into account. The use of effective vertices [221[23] ensures the exact gauge invariance of 
calculated amplitudes despite the off-shell initial partons. 

The off-shell amplitude of subprocess (1) reads: 

Ml = eege>^ip)vs,ip 2 )'['HqqUs 2 {Pi) 7-^-77-“T- ^^ ^ 

iPi +P2)-rnjj- ituhTh 

V 2 ^ F (^ 2 ^ 2 )Fq*q*^'Ur '2 (^^l^l), 

s — — tmzi- z 


X 


ZZH 




iki k2)'^iki -|- ^ 2 )^ 


nir 


where e and Cq are the electron and incoming quark (fractional) electric charges, is the 
polarization 4-vector of produced Z boson, s = (/ci -|- /c 2 )^, ki = XiU -|- kix (with i = 1 
or 2), li and I 2 are the 4-momenta of colliding protons, Xi and X 2 are the corresponding 
momentum fractions and T// is the full decay width of Higgs boson, mz and T^ are the 
mass and full decay width of Z boson. We will take the propagators of intermediate 
Z and Higgs bosons in the Breit-Wigner form to avoid any artificial singularities in the 
numerical calculations. The fermion and gauge boson to Higgs vertices are usual: 


T 


Hqq - 


e niqi 
sin 26w mz ’ 


(5) 


-pliu 

^ ZZH 


■g^'^mz, 


sin 20 ^ 1 / 

where rriqi is the mass of produced quark or antiquark, 9\y is the Weinberg mixing angle. 
We will neglect the masses of initial quarks compared to the masses of hnal state par¬ 
ticles but keep their non-zero transverse momenta: 




— k^rp ^ 0 , 


^ 2 T 


— k\rj^ 7 ^ 0 . 


The effective vertex F^*-*^ which describes the effective coupling of off-shell quark and 
antiquark to Z boson reads |22l[23] (see also 


^ q*q*Z 


Y - ki 




-h 


1^2 


ik-k2) \l2-ki) 


(C« - cp''), 


(7) 


where Cy and C\ are the corresponding vector and axial coupling constants. The effective 
vertex T^*q*z satisfy the Ward identity Fg*-*^(fci -|- /c 2 )^ = 0. The off-shell amplitude of 
subprocess (2) reads: 


M 2 


= eCqigYsY^e^ip) Vs^ip2)FYus2iPi) 
+ eeg^V(5“¥e^(p) VskP2)F^^Us^ipi) 


iki + /C 2 )^ 


iPi + ^2)^ 


vYx2l2)Y^*q*gUr2iXlll) + 
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( 8 ) 
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where e^/ is the produced quark (fractional) electric charge, g is the strong charge, a and 
b are the eight-fold color indexes, and 


Fr 



p + pi + rug^ ^ ^ -p-p2 + 

(p-Fpi)2 - {-p - P2y - 


(9) 
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Fr = r 
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p(+) M 
^ F<iZ 


{k2,p)- 


-k2 -p 


ri*™ +P 2 ) + A^^(/ci, -k 2 ,p,Pi+P 2 ). 


( 10 ) 


{-k2-py 

The on-shell quark coupling to the Z boson is taken in a standard form: 


- cw). 


The effective vertices can be written as 


/A* 

pM _ _ t 1 _2 

r q*q*g — 7 /i 7 \ ^2 , 


(/l-fcs) "(/2-fcl)' 


( 11 ) 


( 12 ) 


= (13) 

= (14) 

The corresponding couplings of the off-shell quark or antiquark to usual on-shell quark 
and Z boson are constructed as it was done earlier I2S1: 


r't’K*:. 9 ) = r™")*;. (15) 

The induced term A^^(/ci, k 2 ,qi,q 2 ) has the form [27] : 


A>^'^{ki,k2,qi,q2) 


ki 


‘K 


ill ■ h)iq2 ■ h 


-|- k^ 


l^l^2 


{qi ■ l2)iq2 ■ k) 


(16) 


The summation on the produced Z boson polarizations is carried out with the usual 
covariant formula: 

Ei''(p)i“'(p) = -9'‘" + S^- (17) 

'"'Z 

In according to the fc-r-factorization prescription [IZliTS], the summation over the polar¬ 
izations of incoming off-shell gluons is carried with ^ e^e*^ = k^k^/k|,. In the collinear 
limit, when |k(r| —)■ 0, this expression converges to the ordinary one after averaging on the 
azimuthal angle. In according to the using of the effective vertices, the spin density matrix 
for off-shell spinors in initial state is taken in the usual form Y2'^{^ih)u{xili) = xpi + m 
(where i = 1 or 2 and we omited the spinor indices). Further calculations are straightfor¬ 
ward and in other respects follow the standard QCD Feynman rules. The evaluation of 
traces was performed using the algebraic manipulation system FORM [2S] . We do not list 
here the obtained lengthy expressions because of lack of space. The off-shell amplitude of 
gluon-gluon fusion subprocess (3) was derived in our previous paper [29] (see also [30]). 

The cross section of any process in the ^(r-factorization approach is calculated as a 
convolution of the off-shell partonic cross section and the unintegrated, or transverse 
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momentum dependent (TMD), parton densities in a proton. The cross sections of sub¬ 
processes (1) and (2) read: 


a = 


Y. 


2567r3(xiX2s)2 


\M 


1 , 2 | 


X 


X/g(xi, h- )d\<.ijnd\i.2jndpirpP2rpdydyidy2 


d(j)i d(f)2 d'lpi dtl)2 
2tt 2tt 2ti 271 


(18) 


where fq(xi, is the TMD quark density in a proton, y is the rapidity of produced 

Z boson, s is the total energy, pi^, p 2 T, ?/i, y 2 -i V'l cind -02 are the transverse momenta, 
rapidities and azimuthal angles of final state quarks, respectively. The incoming quarks 
have azimuthal angles 0i and 02- The cross sections of subprocess (3) can be written as: 


a 


2567r3(a;ia;2s)' 


■\M. 




(19) 


where /^(xi, k|p,is the TMD gluon density in a proton, and is the off-shell am¬ 
plitude of subprocess (3). 

Concerning the TMD parton densities in a proton, we concentrate on the approach 
based on the CCFM evolution equation [3T]. The CCFM parton shower, based on the prin¬ 
ciple of color coherence, describes only the emission of gluons, while real quark emissions 
are left aside. It implies that the CCFM equation describes only the distinct evolution of 
TMD gluon and valence quarks, while the non-diagonal transitions between quarks and 
gluons are absent. Below we use the TMD gluon and valence quark distributions which 
were obtained [521E3] from the numerical solutions of the CCFM equation (namely, set 
AO). Following to [3l], we calculate the TMD sea quark density with the approximation, 
where the sea quarks occur in the last gluon-to-quark splitting. At the next-to-leading 
logarithmic accuracy Q(s(q(s Inx)", the TMD sea quark distribution can be written as 
follows [51]: 


T(sea) 


(x,k|,/i 2 ) 



X 


qt. q^, y), 


( 20 ) 


where 2 ; is the fraction of the gluon light cone momentum carried out by the quark, and 
A = kr — zqjn. The sea quark evolution is driven by the off-shell gluon-to-quark splitting 
function Pgg( 2 ;, q|,, A^) [35] : 




R 


A^ -I- 2;(1 — z)q^^ 


(1 - zf -1- 2 ;^ -1- 4^^(1 - z) 


2 


( 21 ) 


where Tr = 1/2. The splitting function Pqg{z, q^, A^) has been obtained by generalizing 
to hnite transverse momenta, in the high-energy region, the two-particle irreducible kernel 
expansion (SB]. It takes into account the small-x enhanced transverse momentum depen¬ 
dence up to all orders in the strong coupling constant, and reduces to the conventional 
splitting function at lowest order for |q(r| —t 0. The scale is defined 071 from the an¬ 
gular ordering condition which is natural from the point of view of the CCFM evolution: 
/i2 = AV(l-z)2 + q^/(l-z). 

Other essential parameters were taken as follows: renormalization scale /i^ = rn|-|-py, 
factorization scale = s + Q,^ (with Qt being the transverse momentum of initial 
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parton pair), beauty quark mass mb = 4.75 GeV, mz = 91.1876 GeV, mn = 125 GeV, 
Tz = 2.4952 GeV, Th = 4.3 MeV, sin^ 9w = 0.23122 and we use the LO formula for the 
strong coupling constant with Uf = 4 active quark flavors at Aqcd = 200 MeV, so 

that as{m?z) = 0.1232. To take into account the non-logarithmic loop corrections to the 
production cross sections, we apply the effective 77-factor, as it was done in [381139] : 


K = exp 



a,(/i2) 2 

-TT 

27r 


( 22 ) 


where color factor Gp = 4/3. A particular scale choice = pj, m/ was proposed [3811^ 
to eliminate sub-leading logarithmic terms. Note we choose this scale to evaluate the 
strong coupling constant in (22) only. Everywhere the multidimensional integration have 
been performed by the means of Monte Garlo technique, using the routine VEGAS [40] . 
The corresponding G-|--|- code is available from the authors on requeslU. 

We now are in a position to present our numerical predictions. The differential cross 
sections of associated Zbb production in pp collisions as a function of M, the invariant 
mass of final beauty quarks, and Z boson transverse momentum at y/s = 8 and 14 TeV 
are shown in Figs. 2 and 3. The solid, dashed and dash-dotted histograms correspond to 
the contributions from the subprocesses (1), (2) and (3), respectively. There is no any cuts 
applied. One can see that the associated Higgs (decaying into the 66 pair) and Z boson 
production cross section lies below the QGD backgrounds by several orders of magnitude 
in a whole pt range, but peaks near Higgs mass. To increase the relative contribution 
from Higgs signal, we repeated the calculations in the restricted region of M, namely 
120 < M < 130 GeV (see Fig. 4). We found that here the associated Higgs and Z boson 
production gives a sizeble contribution to the Zbb cross section at high Z boson transverse 
momenta. So, at y/s = 8 TeV it practically coincides with the leading contribution from 
the gluon-gluon fusion subprocess at pt > 200 GeV. At y/s = 14 TeV, it lies below the 
latter. However, these contributions are almost comparable at pt > 300 GeV. With the 
expected LHG luminosity of about 40 fb“^, our estimation gives 400 — 500 events (with 
beauty quarks originating from the Higgs boson decays) for both energies, 8 and 14 TeV. 
Therefore, the possibility for the experimental detection of Higgs signal appears in the 
kinematical region dehned above. 

A special opportunity to detect the decays of scalar Higgs bosons can be provided 
by the investigations of different angular correlations between the hnal state particles. 
As an example, we calculated the distributions on the angle 6 between the produced Z 
boson and 6 -quark in the Gollins-Soper frame (where z axis is defined with respect to the 
bisector of colliding protons in the bb rest frame), and on the azimuthal angle difference 
A(j) between the hnal beauty quarks in the pp center-of-mass frame. The results of our 
calculations performed near Higgs boson peak (with 120 < M < 130 GeV) are shown in 
Figs. 5 and 6, where an additional cut pr > 200(300) GeV is applied at y/s = 8(14) TeV. 
As it was expected, the isotropic decay of scalar Higgs particle H ^ bb greatly differs from 
the angular distributions predicted by the off-shell amplitudes of subprocesses (2) and (3). 
Moreover, the beauty quarks, originating from the Higgs boson decay, populate mostly at 
low A(j) (see Fig. 6 ), whereas the leading QGD background, as given by the gluon-gluon 
fusion subprocess (3), has more hat A(j) distribution. So, the diherent angular correlations 
between the hnal state particles in the associated Zbb production are very sensitive to the 
source of bb pairs, and therefore future experimental investigations of such observables at 
the LHG with increased luminosity can give a clear information about Higgs signal. 

Finally, we study the size of theoretical uncertainties of our calculations connected 
with the hard scale. As usual, in order to estimate these uncertainties we vary the scales 

^lipatov@theory.sinp.msu.ru 
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by a factor of 2 around their default values. Also, we use the CCFM set A0+ and AO— 
instead of the default TMD gluon density AO. These two PDF sets represent a variation 
of the hard scale involved in (18) and (19). The A0+ stands for a variation of 2/i, while 
set AO— reflects /r/2. We observe a deviation of about 50% with both A0+ and AO— 
sets (see Fig. 7) for the QCD background (as given by the sum of gluon-gluon fusion 
and quark-antiquark annihilation subprocesses considered above). Despite the relatively 
large band of uncertainties, the latter does not change our conclusions. Additionally, to 
investigate the role of higher-order QCD corrections, in Fig. 7 we presented the results for 
the QCD background obtained in the framework of collinear QCD factorization at LO. 
We hnd that in the kinematical region where the possible Higgs signal could be observed 
these corrections are important. 

To conclude, in the present note we applied the fc^-factorization approach of QCD to 
study the possibility to detect the scalar Higgs boson decay if —)■ 66 in the associated Z 
and bb production at the LHC. Our consideration was based on the off-shell production 
amplitudes of q*q* —)■ ZH —)■ Zq'q', q*q* —)■ Zq'q' and g*g* —)■ Zq'q' partonic subprocesses 
supplemented with the CCFM dynamics of parton densities in a proton. The main part 
of higher-order QCD corrections (corresponding to the log 1/x enhanced terms in pertur¬ 
bative series) is effectively taken into account in our consideration. We demonstrated that 
the if —)■ 66 signal can be observed at large transverse momenta near Higgs boson peak 
despite the overwhelming QCD background, and pointed out an important role of angular 
correlations between the produced Z boson and 6-quarks. The gauge invariant off-shell 
amplitudes of q*q* —)■ ZH — Zq'q' and q*q* —)■ Zq'q' partonic subprocesses, calculated 
for the hrst time, can be implemented in a different Monte Carlo event generators, like 
as, for example, cascade [TT] . 
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c) d) 

Figure 1: Examples of Feynman diagrams corresponding to —)■ ZH —)• Zbb (a), 

g*g* —)■ Zbb (b) and q*q* —?■ Zbb (c, d) subprocesses. Full set of diagrams can be obtained 
by permutations of quark, gluon, Higgs and Z boson lines. 




Figure 2: The associated Z + bb cross sections in pp collisions calculated as a function of 
invariant mass of bb quarks at ^/s = 8 TeV (left panel) and ^/s = 14 TeV (right panel). 
The solid, dashed and dash-dotted histograms correspond to the contributions from the 
q*q* —)■ ZH —)■ Zbb, q*q* —)■ Zbb and g*g* —)■ Zbb subprocesses, respectively. No cuts is 
applied. 
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Figure 3: The associated Z + bb cross sections in pp collisions calculated as a function 
of Z boson transverse momentum at ^/s = 8 TeV (left panel) and y/s = 14 TeV (right 
panel). Notation of all histograms is the same as in Fig. 2. No cuts is applied. 




Figure 4: The associated Z + bb cross sections in pp collisions calculated as a function 
of Z boson transverse momentum at ^/s = 8 TeV (left panel) and y/s = 14 TeV (right 
panel) at 120 < M < 130 GeV. Notation of all histograms is the same as in Fig. 2. 




Figure 5: The associated Z + bb cross sections in pp collisions calculated as a function 
of angle 6 between the produced Z boson and beauty quark in the Collins-Soper frame 
at y/s = 8 TeV (left panel) and y/s = 14 TeV (right panel) at 120 < M < 130 GeV. 
An additional cut pt > 200(300) GeV is applied for yfs = 8(14) TeV. Notation of all 
histograms is the same as in Fig. 2. 
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Figure 6: The associated Z + bb cross sections in pp collisions calculated as a function of 
azimuthal angle difference A(j) between the produced beauty quarks in the pp center-of- 
mass frame at \/s = 8 TeV (left panel) and \fs = 14 TeV (right panel) at 120 < M < 
130 GeV. An additional cut pt > 200(300) GeV is applied for ^/s = 8(14) TeV. Notation 
of all histograms is the same as in Fig. 2. 
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Act) [rad] A4) [rad] 

Figure 7: The associated Z + bb cross sections in pp collisions calculated as a func¬ 
tion of Z boson transverse momentum pt-, angle 9 and azimuthal angle difference A0 at 
^/s = 8 TeV (left panel) and y/s = 14 TeV (right panel) at 120 < M < 130 GeV. The solid 
and dash-dotted histograms correspond to the QCD background (sum of the gluon-gluon 
fusion and quark-antiquark annihilation subprocesses) calculated in the framework of 
factorization approach and collinear approximation of QCD at LO, respectively. The up¬ 
per and lower dashed histograms correspond to the scale variations in the fcr-factorization 
predictions, as it is described in the text. The dotted histograms correspond to the con¬ 
tributions from the q*q* —)■ ZH Zbb subprocess. An additional cut pt > 200(300) GeV 
is applied for ^/s = 8(14) TeV in the 6 and A0 distributions. 


13 




















